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Abstract. We show that the spectrum of the Schrodinger operator on a finite, metric 
graph determines uniquely the connectivity matrix and the bond lengths, provided that 
the lengths are non-commensurate and the connectivity is simple (no parallel bonds 
between vertices and no loops connecting a vertex to itself). That is, one can hear the 
shape of the graph! We also consider a related inversion problem: A compact graph 
can be converted into a scattering system by attaching to its vertices leads to infinity. 
We show that the scattering phase determines uniquely the compact part of the graph, 
under similar conditions as above. 
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1. Background and notations 

The question "Can one hear the shape of a drum?", was posed by Marc Kac ]I] as 
a paradigm example for a class of problems which is of fundamental importance in 
many physical applications: Given the quantum spectrum, can one deduce uniquely 
the basic interactions or the geometric constraints which specify the system? In Kac's 
original paper, this inversion problem is formulated for Laplacians on compact domains 
with boundary conditions (billiards), and in spite of approximately fifty years of active 
research, the complete answer is still an enigma. If the billiard boundary is allowed to 
have corners, one can draw different yet isospectral billiards [^|, [§. Isospectral microwave 
cavities were also constructed H using the same models. Recently, it was shown that 
Kac's question is answered in the affirmative for simply connected doamins with analytic 
boundaries with some symmetry restrictions || . For boundaries in intermediate 
classes of smoothness the answer is not known. The existence of isospectral systems 
was investigated for Laplacians on closed Riemannian manifolds || and for discrete 
Laplacians which are formed by the connectivity matrices of graphs . In both cases, 
elaborate techniques were devised to identify large sets of different, but isospectral 
systems. However, if the domains are analytic surfaces of revolution, the spectrum 
determines the manifold uniquely ||. 

In the present note we address yet another class of problems which go under the 
name "quantum graphs". The Schrodinger operator consists of the one- dimensional 
Laplacian on the bonds, endowed with boundary conditions on the vertices, which 
ensure that the spectra are real and discrete, (see P for a detailed discussion and 
list of references). These systems recently attracted the attention of the quantum 
chaos community, since, in spite of their structural simplicity, their spectra display 
the complexity and statistics which characterize generic quantum chaotic systems such 



as billiards. Moreover, an exact trace formula |Tl], [L2| can be written down in terms 
of the periodic orbits (PO) on the graph, which is analogous to the semiclassical trace 
formula for chaotic systems. So far, this exact trace formula paved the way to the study 
of spectral statistics by applying combinatorial methods to periodic orbit expansions 



|13| , |14}| . Here we shall show that it also provides the key to the affirmative answer to 
Kac's question, when applied to graphs. 

A graph Q consist of V vertices connected by B bonds. The V x V connectivity 
matrix is defined by: 

Cjj = number of bonds connecting the vertices i and j . (1) 

A graph is simple when for all i, j : G [0, 1] (no parallel connections) and = 
(no loops). The valence of a vertex is Vi = Y^j=i^i,j an d the number of bonds is 



B = \ j = iCij. We denote the bonds connecting the vertices i and j by b = 
The notation [i,j] and the letter b will be used whenever we refer to a bond without 
specifying a direction. Hence, b = [i,j] = [j,i]. To any vertex i we can assign in a 
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unique way the set of bonds which emanate from it: 

s {l) = {all bonds [i, k] : C itk = 1} ; card[s^} = v t . (2) 

We shall refer to sW as a topological star, since its definition does not require any metric 
information. For simple graphs, 

C id = C jti = card[s {i) n s {j) ] . (3) 

Directed bonds will be denoted by d — and we use the convention that the bond 

is directed from the first to the second index. The round brackets, and the letter d 
are reserved to denote the directed bonds. The reverse bond to d — is denoted 
by d = (J,i). In analogy to (||), we define the sets of directed bonds outgoing from or 
incoming to i: 

s (h+) _ jgji di rec ted bonds (i, k) : = 1} 

S ( J ■>-) = {all directed bonds (k, i) : C k ,i = 1} , (4) 

respectively. Note: card[s^' + ^] = card[s^'^} = i>j. A directed bond d' is defined to 
follow the directed bond d at the vertex i if d G s^~^ and d' G s^' + \ It is convenient 
to define a directed connectivity matrix F, with = 1 if d follows d', and F djd i = 
otherwise. Any directed bond d = (i,j) belongs to the two sets, s^ ,+ ) and s^~\ and 
C id = card[s^ n s^'-)] = card[s^> + ^ n s^' - )]. 

We shall now define the Schrodinger operator on Q, and collect a few facts which 
will be used in the sequel. For a detailed exposition see e.g., 0. We assign the natural 
metric to the bonds, and each bond is endowed with a length The Schrodinger 
operator consists of 



• d / 2 
-l-T- - A h 

ax 



(5) 



on all the bonds, and At, are constants which are introduced to break the symmetry 
between the counter propagating waves on the bond. The differential operator is 
supplemented by boundary conditions which are imposed in the following way. The 
most general solution of the Schrodinger equation on any bond b = [i,j] is 

ip b {x) = a d e i{k+Ab)x + a d e- i{k - Ab)x ; V b G s {i) , (6) 

where d, and d denote the two directions on the bond b, and a d , a d are yet unspecified 
complex numbers. The boundary conditions impose the relations 

W G : a d = °% a * > ( 7 ) 

where the vertex scattering matrix is a unitary, V{ x V{ matrix, can be either an 
arbitrary k independent unitary matrix, or of the general class of (possibly k dependent) 
matrices derived from matching conditions at the vertices ||, [TD|. The matrix crW will 
be called properly connecting if none of its matrix elements vanishes. The transition 
probabilities W dd , = \o~ d % \, | 2 appear in the classical dynamics analogue, which is a 
random walk model on the directed bonds, with transition probabilities W^,. As 
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an example, consider the most commonly used and discussed quantum graphs, where 
the vertex scattering matrices are derived by imposing Neumann boundary conditions 
at the vertices. In this case a d \, = —5 d , j + ^ , and but for the trivial vertices with 
Vi = 2, these scattering matrices are properly connecting. 

Combining the boundary conditions for all the vertices, results in a secular equation 
for the wavenumbers k n which consists of the spectrum of the Schrodinger operator. The 
secular equation reads, 

det(/-5(jfe)) =0 (8) 

where S(k) is a IB x IB unitary matrix, whose rows and columns are labeled by the 
directed bond labels. It is defined as 

S d ,Ak) = F d)d ^ k+A ^a d % , (9) 

The index i denotes the vertex at which d follows df, and L d = L d , is the bond 
length which is independent of the direction of propagation. However A d = —A d , 
which distinguishes between the directions of propagation. S(k) can be interpreted as 
a quantum evolution operator describing the scattering of waves with wave number k 
between connected bonds. The wave which scatters at the vertex i from d' G s %, ~ to 
d G s l,+ with an amplitude <r d \,, gains the phase (k + A d )L d during the propagation 
along the outgoing bond d. With the next application of S(k) the wave scatters again 
on the vertex to which d is directed. If all a® matrices are properly connecting, the 
waves on Q propagate between all the topologically connected bonds. The unitarity 
of the matrices guarantees that S(k) is unitary for real k. This is also why the 
Schrodinger operator is self-adjoint and its spectrum is pure-point and unbounded on 
the real line. 

The secular equation (|j), with the form (|9D for S(k) leads naturally to the exact 
trace formula |Tl], [12[ 

d(k)= J2 5 ( k ~ k n) 
n 

= - + ^- Y (A p e ikl ? + A*e- ikl -) (10) 

P 

where £ = J2b=i Lb is the total length of Q. The sum extends over all the PO's on 
the graph, and goes over primitive PO's as well as their repetitions. The length of a 
PO is denoted by l p . The coefficient A p is a product over all the scattering amplitudes 
<J dd i, encountered along the PO, times the length of the primitive PO of which p is a 
repetition. It is also endowed with a phase factor exp(z A d L d ). The length of a PO is 

B 

1 p = J2 i ( b p) Lb > * G ^ ( n ) 

6=1 

The length spectrum is the Fourier transform of ( |I0"D . 

n p 
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2. Hearing the shape of the graph 

After reviewing the necessary background, we are in position to formulate the answer 
to Kac's question for graphs. 

Theorem 1: The spectrum of a Schrodinger operator on a metric graph determines 
uniquely the graph connectivity and the length of the bonds, provided that 

- the graph is finite and simple, 

- the bond lengths are rationally independent, 

- the vertex scattering matrices are properly connecting. 

Before presenting the formal proof, we shall sketch its main idea: The length 
spectrum(|l^) is constructed from the Schrodinger spectrum. Its singularity at I = 
provides the total length of the graph £, which, in turn, gives an upper bound to the 
individual bond lengths. The conditions that the lengths are rationally independent, 
and that whenever d follows d! at vertex z, a^ d , ^ 0, ensure that the lengths of all the 
PO's which are consistent with the graph connectivity, appear as 5 singularities in the 
length spectrum, and their corresponding coefficients A p do not vanish. Moreover, PO's 
p, p' which traverse the bonds different numbers of times ( ^ ' for some b in 
(PI)) have distinct lengths. We use these facts to isolate two classes of periodic orbits. 
The 2-PO's of the type [d, d] provide the lengths of the bonds. The 4-PO's of the type 
[d®, dP, d^ m \ d^\ for a fixed / and for all possible m are used to identify the topological 
star at the vertex i where d^ follows d"'. Once all the topological stars s^' are found, 
we apply (|3]) to get the connectivity. Thus, the shortest, and structurally simplest PO's 
are used to extract the parameters which determine the "shape" of Q. 

The proof of the theorem proceeds as follows: 

i) . Use the spectrum to construct the sum dl2|) . Its singular support R consists of I = 
and the infinitely many points l p of the form (|TT|) on the real axis. The weight of the 
singularity at I — is 2£. 

ii) Generate a finite set P which consists of all the strictly positive lengths l p G R, which 
are less than 2C. Exclude any length which can be expressed as an integer multiples 
of any other length, with a multiplier larger than 1. P is the (finite) set of lengths of 
primitive PO's on Q, with lengths in (0,2/2]. 

Hi) Define the basis Q C P, which consists of the minimum number of lengths X q G P 
such that 

card[Q] 

\fl p E P , l p = - n g X q , n q eN (13) 

9=1 

P is a finite set, and therefore Q can be constructed in a finite number of steps. Because 
of the rational independence of the bond lengths, the basis is composed of the lengths 
of 2-PO's [d,d], with Aq = 2L q where L q are the lengths of the bonds. card[Q] = B is 
the number of bonds, and X q = 2£. 

iv) In what follows we shall construct metric stars which will be shown to be in one-to- 
one correspondence with the topological stars (@): Define two lengths \ q and \ q > G Q to 
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be connected if X q + X q i G P. A length in Q cannot be connected to itself since repetitions 
where excluded from P. The lengths in Q which are not connected to any other length, 
correspond to disconnected bonds, which are identified this way, and excluded from 
the susequent analysis. Given a length A r G Q, consider all the clusters C m (X r ) C Q 
satisfying the following requirements: (a) A r G C m (X r ). (b) All the lengths in C m (A r ) 
are pairwise connected, (c) C m (X r ) is maximal - adding another X s G Q, will violate 
(b). Denote the number of C m (A r ) by N(X r ). Metric stars are special clusters defined as 
follows. If for some X s , N(X S ) = 1 we shall define {A s }, as a metric star consisting of a 
single element. If card[C m (X s )] = 2 then it is a metric star with two elements. All other 
C m (X r ) are metric stars if for any three elements A;, X m , X n G C m (A r ), Xi + X m + X n G P. 
P is pruned of repetitions. Hence, this condition excludes repeated triangular PO's from 
being considered as stars. Using all the A r G Q as reference lengths, we construct all 
the different metric stars in Q. The bond lengths are rationally independent, thus, the 
metric stars S^&re in 1 <-> 1 correspondence with the topological stars (|2]). The proof 
is completed by writing the connectivity matrix in analogy with (^): 

C itj = card[S® n S®] ; L M = ^ (S® n S^) . (14) 



The length spectrum can be used to get more information about the graph |15|]. The 
coefficients of the 5 singularities provide constraints on the values of the elements 
of the vertex scattering matrices. Assuming that the a® are symmetric, one can 
determine the transition probabilities W d \, from the coefficients of the lengths of three 
classes of simple PO's. Denote b = = [i, k],b" = [i,l] and correspondingly 

d = (i,j),d' = (i,k),d" = They belong to the same star. The three types of 

PO's are [d,d\, [d',d',d,d], and [d" , d" , d' , d' , d,d], with lengths 2L b , 2{L b + Ly), and 
2(L b + Ly + Lyr), with corresponding coefficients 

\a = 2 L » "fft < 15 > 

A lti , dA = 2(L b + L i ,)af^fy<^ (16) 
•^[d",d",d',d',d,d] = 4 ( L b + L v + Lb") x (17) 

(i) (*0 W (?) (*) (0 
d>'d d,d d,d' d'd' d'd" d" ,d" 



Dividing (|1^) and (|T7| ) by ( |T5| ) we get ratios of the matrix elements of a single crW of 
the form 

A W (*) 

•A-[d'd>,d,d) 2L b L b > a d',d a d,d' / 1 x 



Ratios of products of three matrix elements can also be formed. Using this information 
one can construct the matrices uniquely, up to a phase factor since the ratios of 
the type (|18D are invariant under a right and left multiplication by arbitrary diagonal 
matrices. Hence one obtains only absolute values, in other words, the transition 
probabilities W& d ,. When <j® are real, L b A b m.od2'K can also be computed from the 
weights of simple PO's. 
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3. Application to non-compact (scattering) graphs 

A related inversion problem, is encountered when the compact graph is turned into a 



scattering system by coupling to its vertices leads (bonds) to infinity |jIE| . It is assumed 
throughout that each vertex can be attached to at most one lead. The definition of 
the Schrodinger operator can be extended naturally to the non-compact case, and the 
coupling to the leads is achieved by modifying the vertex scattering matrices to include 
coupling to the leads. Let N be the number of vertices to which leads are attached, and 
the leads are denoted by the index of the vertices to which they are connected. Then, 
for every wavenumber k one can find a solution to the Schrodinger equation which on 
the leads takes the form 

iPiix) = h{k)e- ikx + O^k)^ (19) 

and the outgoing amplitudes Oi(k) are related to the incoming amplitudes h{k) by 

N 

Oi{k) = Y,TiAk)h{k) (20) 
v=x 

where T(k) is the NxN scattering matrix. The conservation of flux ensures the unitarity 



of T(k). An explicit expression for the scattering matrix is given in [I6|. Here, we shall 
be concerned with the scattering phase defined by 

e(fc) = -^log[detT(fc)] (21) 

ZTCt 

The scattering phase is the counterpart of the spectral counting function for compact 
systems. This can be best seen by considering the function dn(k) = d ^1 ■ ft assigns 
a normalized Lorentzian to each scattering resonances, and it is the counterpart of the 
spectral density fllPl). An exact trace formula for dn(k) was derived in [pij. It reads, 

In 



7Z {k-k n f + 1 i 
i 

7T ' 2"7T 



L — y"(A p e lkl p + A* v e- ikl i>) (22) 



V 

where k n — vy n are the poles (resonances) of T(k) in the complex k plane. As in the 
compact case, £ = J^6=i ^ b ^ s ^ e total length of the bonds in the compact part of the 
graph. The sum extends over all the PO's (periodic orbits) on the compact part of the 
graph, and goes over primitive PO's as well as their repetitions. The length of a PO is 
denoted by l p . The amplitude A p is a product over all the scattering amplitudes <5^,, 
encountered along the PO, times the length of the primitive PO of which p is a repetition. 
It is also endowed with a phase factor ^ A^L^. The vertex scattering matrices a^' in 
the space of the directed compact bonds are sub-unitary (| detail < 1). To complete 
0"^ to a unitary matrix a row and a colomn should be added, to take into account the 
coupling of the compact bonds to the lead. For this reason, the rows of the classical 
transition matrix Wj d , = |5"^/| 2 do not sum up to unity, which take into account the 
flux which escapes through the leads. The same mechanism is responssible also to the 
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Figure 1. Two isospectral graphs. The bond lengths are expressed in terms of the 
two arbitrary lengths a and b. 

reduction of \A P \ relative to the corresponding \A P \. This accelerates the convergence 
of ( ^2|) and pushes its poles away from the real axis. A length spectrum similar to (|l|) 
is obtained from the Fourier transform of (p2[), and its singular support consists of 
and the length spectrum of the compact part of the graph. The following result can be 
proved using the same argument as above: 

Theorem 2: The scattering phase Q(k) of a non-compact metric graph, determines 
uniquely the connectivity and the bond lengths in the compact part of the graph, 
provided that 

- the graph is composed of a finite and simple compact part, with leads to infinity. A 
vertex can be coupled to at most one lead. 

- the bond lengths in the compact part are rationally 
independent, 

- the vertex scattering matrices are properly connecting. 

The proof follows verbatim the proof of theorem 1. 

4. An example of isospectral graphs 

We have shown that the spectral density (for compact graphs) and the resonnance 
density (for scattering graphs) determine uniquely the bond lengths and the the 
connectivity, provided that the bond lengths are rationally independent, and the vertex 
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scattering matrices are properly connecting. We shall now construct a non trivial 
example of two different, yet isospectral graphs, obtained when the requirement of 
rational independence is relaxed. 

Extending the geometric construction of || we produced the pair of graphs shown 
in Figure 1. Their bond lengths are rational combinations of two arbitrary lengths a and 
6, as shown in the figure. Neumann boundary conditions are imposed at the vertices. 
The spectra of the two graphs are the zeros of the corresponding secular functions (f|), 
which for the present cases read: 

Z {I) (k) =tg(2 (a + b) k) + (23) 

2 tg(a k) + 2 tg(6 k) + tg((2 a + 6) k) + tg((q + 2b) k) 
1 - (2 tg(a k) + tg(6 k)) (tg(6 k) + tg((2 a + b) k) + tg((a + 2 6)*)) 

Z (II) (k)=tg(2ak) (24) 

2 tg(a k) + 2 tg(6 k) + tg((a + 2 b) k) + tg((2 a + 3 b) k) 
1 - (tg(a fc) + tg(6 k) + tg((a + 2 6) k)) (tg(o fc) + tg(6 k) + tg((2 a + 3 b) k)) 
It is not an easy matter to show directly that these two different functions vanish at 
exactly the same values of k. However, the theorem which underlies their construction 
|| guarantees this fact. 
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